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ADDITIVITY  OF  RETINAL  DAMAGE  FOR  MULTIPLE-PULSE 
LASER  EXPOSURES 


INTRODUCTION 
MilUary  Relevancy 

The  ever-increasing  use  rf  lasers  In  military,  medical,  and  industrial 
environments  presages  increasing  incidents  of  personnel  exposed  to  harmful 
levels  of  laser  radiation.  While  a  broad  data  base,  accumulated  over  the 
past  two  decades,  defines  ocular  damage  thresholds  for  single-pulse  expo¬ 
sures  from  either  continuous-wave  (cw)  or  pulsed  sources,  relatively  few 
studies  have  examined  the  effects  of  multiple-pulse  or  repeated  exposures. 
However,  a  field  scenario  can  easily  be  envisioned  in  which  the  dynamics  of 
laser  beams  and  targets,  both  in  motion,  dictate  that  a  sequence  of  expo¬ 
sures  is  as  likely  as--or  more  likely  than--a  single-exposure  episode.  If  a 
repetitively  pulsed  laser  is  Involved  with  a  pulse-repetition  frequency 
greater  than  10  Hz  (up  to  megahertz),  a  few  or  m?ny  pulses  may  strike  the 
retina  with  each  exposure  encounter.  Because  of  the  task  or  mission  at 
hand,  the  targeted  or  accidentally  exposed  personnel  could  have  a  tendency 
(if  not  a  compelling  reason)  to  look  in  the  direction  from  which  the  laser 
beam  emanates,  so  that  areas  of  exposed  retinal  tissue  could  be  concentrated 
in  the  macular  area.  Under  such  a  scenario,  the  probability  of  a  given 
patch  of  retinal  tissue  receiving  multiple  exposures  might  be  considerable, 
and  would  Increase  with  the  time  that  personnel  remained  vulnerable  to  laser 
exposures.  Thus,  in  evaluating  potential  ocular  effects  of  lasers, 
achieving  an  understanding  of  the  cumulative  effects  of  repeated  exposures 
is  of  significance. 


Background 

The  genesis  of  thermal  damage  additivity  studies  is  the  work  of 
Henriques  (1)  and  Moritz  and  Henriques  (2).  In  their  studies  of  thermal 
injury  to  skin,  they  observed  additivity  of  damage  from  exposures  over 
periods  of  several  hours;  and  they  defined  the  Henriques  damage  integral  as 
a  measure  of  extent  of  damage  as  a  f”nctioii  of  the  spatial -temporal  tem¬ 
perature  increases  associated  with  a  thermal  in'uilt.  The  Henriques  damage 
integral  has  been  incorporated  into  .hermal  models  of  laser-induced  ocular 
damage  developed  by  Mainster  et  al .  (3)  and  Takata  et  al ,  (4).  The  model  of 
Takata  et  al .  has  been  tested  against  experimental  thresholds  for  ophthal- 
moscopically  visible  lesions  induced  by  single-pulse  laser  exposures,  and 
has  proved  to  be  a  highly  reliable  predictor  of  retinal  damage  threshold  for 
pulsewidths  ranging  from  10“7  to  10^  s  (4).  For  multiple-pulse  exposures, 
the  thermal  model  was  not  tested  In  any  systematic  manner. 
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The  need  for  further  quantitative  studies  of  the  additive  effects  of 
multiple-pulse  exposures  is  perhaps  best  illustrated  by  examining  the  treat¬ 
ment  of  multiple  pulses  in  existing  laser  safety  standards.  The  current  Air 
Force  Occupational  Safety  and  Health  'AFOSH)  standard  (5)  treats  pulses 
separated  by  s  as  totally  independent  events  (i.e,,  zero  additivity). 

The  American  National  Standards  Institute's  (ANSI)  Standard  for  the  Safe  Use 
of  Lasers,  prior  to  a  1986  revision,  required  that  pulses  at  repetition 
rates  of  <1  Hz  be  added  over  a  24-hr  period  (i.e.,  total  additivity)  (6). 

The  revised  ANSI  standard  (7)  defines  multiple-pulse  exposures  only  for 
cases  with  a  pulse  repetition  frequency  (PRF)  >1  Hz  and  thus,  by  omission, 
leaves  exposures  separated  by  H  s  to  be  treated  as  Independent  single-pulse 
events.  While  the  old  ANSI  position  of  total  additivity  over  a  24-nr  period 
may  have  been  unnecessarily  conservative,  the  current  ANSI  and  AFOSH  treat¬ 
ments— with  no  additivity  for  exposures  separated  by  more  than  1  s--may  be 
Inadequate,  in  some  cases,  for  defining  suitable  protection  from  multiple- 
pulse  exposures. 

The  most  recent  revision  of  the  ANSI  laser  safety  standard  suggests  an 
empirically  derived  treatment  for  multiple-pulse  exposures  (7).  According 
to  the  standard,  the  maximum  permissible  exposure  (MPE)  per  pulse  for  a 
train  of  n  identical  pulses  is  equal  to  {n)“l/^  times  the  MPE  for  a  single 
pulse  in  the  train.  This  simple  and  totally  empirical  relationship  provides 
an  acceptable  approximation  to  the  body  of  multiple-pulse  threshold  data 
which  existed  at  the  time  this  relationship  was  first  proposed  (8,9),  Ob¬ 
viously,  one  must  use  great  caution  in  generally  applying  a  relationship 
which  depends  on  only  the  one  exposure  parameter  (n),  and  which  has  not  been 
systematically  tested  (because  of  a  lack  of  experimental  data),  versus  wide 
ranges  of  other  exposure  parameters  (e.g.,  PRF,  pulse  train  length,  beam 
divergence,  etc.). 

The  objective  of  the  current  work  is  to  provide  a  quantitative  basis 
for  evaluating  the  cumulative  nature  of  multiple-pulse  or  repeated  exposures 
for  the  various  exposure  conditions  which  may  be  anticipated  with  existing 
laser  systems.  In  particular,  this  work  is  concerned  with  additivity  of 
repeated  subthreshold  laser  exposures,  where  the  exposure  conditions  ire 
such  that  any  resultant  retinal  lesion  would  arise  from  a  thermal  damage 
mechanism.  Therefore,  this  work  encompasses  very  broad  ranges  of  exposure 
parameters,  but  excludes  long  exposures  (>1  s)  to  short  visible  or  ultra¬ 
violet  wavelengths  where  photochemical  processes  might  result  in  retinal 
damage  at  exposure  levels  Insufficient  to  induce  thermal  damage  (10).  Also 
excluded  are  u’tra-short  pulsewldth  exposures  (of  a  nanosecond  or  less) 
where  photo-acoustic  and  other  non-linear  mechanisms  may  be  dominant  (11). 

Two  basic  tissue  rate  processes  must  be  considered  in  assessing  the 
additivity  of  repeated  pulses.  The  first  is  the  thermal  conductivity  of  the 
tissue.  For  the  retina,  thermal  relaxation  time,  are  on  the  order  of 
milliseconds  (12).  The  second  rate  process  is  that  of  tissue  repair  of 
reversible  damage.  The  existence  of  a  repair  process  is  an  a  priori 
assumption  based  on  observations  that  repeated  low-level  exposures  (e.g., 
ambient  sunlight)  are  not  cumulative  over  indefinite  periods  of  time,  or  we 
would  all  suffer  severe  retinal  degeneration  in  short  order.  Studies  in 
this  laboratory  have  quantitatively  assessed  corneal  and  retinal  repair 
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rates  following  ultraviolet  and  blue-light  exposures  (where  photochemical 
damage  mechanisms  are  involved).  The  time  constants  associated  with  tissue 
repair  were  found  to  be  approximately  2  and  4  days  for  the  cornea  and 
retina,  respectively  (13,14).  That  a  comparable  repair  rate  exists  after 
laser-induced  thermal  insult  to  the  retina  is  uncertain.  However,  the 
organism's  response  may  be  determined  by  the  given  degree  of  tissue  damage, 
without  regard  to  (or  "memory"  of)  the  exposure  parameters  and  specific 
mechanism  which  Induced  that  degree  of  damage. 

Given  the  two  tissue  rate  processes  just  discussed,  considerations  of 
multiple-pulse  laser  exposures  may  relate  to  the  relative  values  of  the  time 
constants  associated  with  these  rate  processes  on  the  one  hand,  and  with  the 
pulsewidth,  PRF,  and  pulse  train  length  on  the  other  hand.  Thus,  for  PRFs 
>>1  Hz,  thermal  relaxation  of  the  expos^^  area  may  not  be  completed  during 
the  interval  between  pulses;  and  ..lie  consequences  of  repeated  pulses  may  be 
additive  in  a  manner  consistent  with  the  higher  and  higher  peak  temperatures 
achieved  with  each  successive  pulse,  A  different  result  may  be  achieved  if 
the  identical  sequence  of  pulses  is  delivered  over  a  longer  time  frame  (with 
PRF  <1  Hz)  where  thermal  relaxation  is  effectively  completed  during  each 
inter'-pulse  interval.  In  such  a  case,  the  peak  temperature  at  any  point 
never  exceeds  that  resulting  from  a  single  pulse  in  the  train. 

Regardless  of  PRF,  if  one  is  concerned  with  evaluating  the  cumulative 
effects  over  periods  of  time  comparable  to  or  longer  than  the  tissue  repair 
rate,  then  the  additivity  of  repeated  exposures  is  diminished  by  the  ongoing 
repair  process.  Oniy  when  the  pulse  train  length  is  much  less  than  the 
tissue  repair  rate  can  the  latter  be  Ignored  in  calculating  the  net  cumu¬ 
lative  effect. 

Ultimately  then,  we  seek  a  quantitative  basis  for  evaluating  tlie 
additive  nature  of  multiple-pulse  or  repeated  exposures  for  the  various 
contingencies  just  mentioned;  i.e.,  for  PRFs  ranging  from  much  less  than  to 
much  greater  than  the  thermal  relaxation  rate,  and  for  pulse  train  lengths 
ranging  from  much  less  than  to  much  greater  than  the  tissue  repair  rate. 


METHOD 

Subjects 

Experimental  subjects  wtn'e  rhtsus  monkeys  (rLcaca  mulatta) .  Potential 
subjects  were  screened  with  a  slit  lamp  and  fup-ins  camera,  and  only  animals 
with  clear  ocular  media  and  normal  fundi  were  used.  In  addition,  no  animal 
was  used  that  had  a  retractive  error  difference  of  more  than  0.50  diopter  in 
any  meridian. 

One  day  before  a  planned  experiment,  2-3  drops  of  atropine  sulfate  (IX 
ophthalmic  solution)  were  introduced  Into  each  conjunctival  sac  to  induce 
cycloplegia.  The  animals  were  tranquili/ed  with  an  Injection  of  ketamine 
hydrochlori  a  (12-14  mg/kg  body  weight,  i.in.),  and  anesthetized  to  effect 
with  sodium  pentobarbital  introduced  through  a  catheter  iri  a  superficial  leg 
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vein.  Throughout  the  experiments,  the  primate  core  temperature  was  moni¬ 
tored  by  means  of  a  rectal  probe,  and  maintained  at  37  by  means  of  a 
thermal  blanket. 


Apparatus 

All  exposures  were  carried  out  with  a  Spectra  Physics  171  krypton  laser 
with  a  cw  output  at  647  nm.  An  electronically  controlled  mechanical  shutter 
and  two  Vincent  Associates  shutter  drive  timers  In  series  provided  a  pro¬ 
grammable  means  of  choosing  the  pulsewldth,  pulse  repetition  rate,  and 
pulse-train  duration.  In  all  cases,  a  lOO-ins  pulsewldth  was  chosen,  but  the 
PRF  and  pulse  train  length  were  varied. 

A  pellicle  beamsplitter  was  positioned  to  deflect  a  fraction  of  the 
shuttered  beam  to  a  silicon  photodiode  detector  and  Photodyne  66XLA  optical 
power-energy  meter.  In  the  power  mode  and  with  the  shutter  open,  this  meter 
monitored  the  cw  output  of  the  laser  before  and  after  each  exposure.  In  the 
energy  mode,  the  meter  measured  the  total  energy  delivered  In  each  shuttered 
pulse  train.  In  either  case,  the  relative  numbers  obtained  from  the  de¬ 
flected  portion  of  the  beam  were  converted  to  absolute  levels.  Incident  at 
the  subject's  cornea,  by  cross-calibrating  against  a  second  photodiode 
detector-Photodyne  meter  combination  where  the  second  detector  was  Inserted 
at  the  position  otherwise  occupied  by  the  subject.  The  cross-calibration 
was  completed  before  and  after  each  exposure  session. 

A  Zeiss  fundus  camera  was  fitted  with  a  swing-out  mirror  which,  when 
out  of  position,  allowed  fundus  camera  viewing  of  the  subject's  retina. 

When  In  position,  the  mirror  deflected  the  Incoming  laser  beam  through  the 
subject's  pupil  and  in  the  direction  colllnear  with  the  optic  axis  of  the 
fundus  camera,  thus  allowing  the  beam  to  be  steered  to  a  pre-selected  spot 
on  the  fundus  according  to  the  positioning  of  the  subject. 

The  krypton  laser  beam  had  a  divergence  of  1  mrad  and  a  diameter  of 
~2  mm  when  incident  at  the  corneal  plane.  To  vary  the  retinal  Image  size  of 
the  beam,  an  ophthalmic  trial  lens  was  placed  In  the  beam  path  at  a  point 
6  in.  from  the  corneal  plane.  For  example,  a  +7  diopter  lens  at  this  posi¬ 
tion  yielded  a  Z00-/im  retinal  image  diameter.  This  was  verified  by  direct 
observation  with  tne  fundus  camera  when  a  pellicle  beamsplitter  replaced  the 
swing-out  mirror  to  permit  fundus  viewing  during  low-level  laser  exposure, 
thus  allowing  the  retinal  image  to  be  measured  by  a  calibrated  reticle 
Incorporated  into  the  fundus  ramera  optics.  When  no  trial  lens  was  used  and 
the  collimated  laser  beam  was  allowed  through  the  pupil,  the  retinal  image 
diameter  was  ~30  jim. 


Laser  Exposures 

Anesthetized  animals  were  mounted  on  an  adjustable  stage  in  front  of 
the  fundus  camera,  and  their  eyes  were  positioned  by  translational  and 
rotational  stage  adjustments.  At  the  beginning  of  each  exposure  session, 
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extraniacular  marker  lesions  were  placed  in  each  eye  to  define  a  3  x  3 
pi.'cular  grid  for  placement  of  experimental  exposures.  Marker  lesions  were 
generated  by  single  100-ras  pulses,  by  using  tne  collimated  647-nm  laser  beam 
with  an  energy  of  8  to  10  mJ.  Marker  lesions,  although  initially  punctate, 
typically  developed  to  ~100-/im  diameter  within  1  hr  postexposure. 

After  the  marker  lesions  developed,  macular  exposures  were  delivered  to 
each  of  the  nine  marker  grid  sites.  The  range  of  doses  was  from  ~0.5  to 
2.0  times  the  estimated  threshold  for  an  ophthalmoscopically  visible  lesion. 
For  the  set  of  exposures  delivered  to  any  given  eye,  only  the  power  of  the 
laser  was  varied,  while  all  other  exposure  parameters  were  held  constant. 
Exposures  were  generally  coinpleted  for  both  eyes  during  a  single  experi¬ 
mental  session.  The  animal  was  maintained  at  an  even  plane  of  anesthesia 
until  the  lesion-no  lesion  ophthalmoscopic  readings  were  completed  at  1  hr 
postexposure. 


Data  Analysis 

The  lesion-no  lesion  data  were  subsequently  analyzed  by  the  method  of 
probit  analysis  (15).  Generally,  data  from  six  eyes  (up  to  54  data  points) 
were  used  for  each  threshold  determination.  A  probit  program  was  used  to 
calculate  the  dose  at  which  there  is  a  50%  probability  of  inducing  damage 
(ED50)  and  the  requested  confidence  limits  (usually  95%)  on  the  ED50.  The 
program  also  generates  a  probability  vs.  dose  curve,  so  that  other  thresh¬ 
olds  (e.g.,  EDio  and  ED90J  are  available  as  desired. 


Thermal  Model  Calculations 


A  retinal  thermal  model  was  used  to  predict  damage  thresholds  for 
comparison  with  the  experimentally  determined  multiple-pulse  thresholds  and 
with  the  empirical  fits  of  multiple-pulse  threshold  data  (8,9).  The  model, 
initially  developed  by  Takata  et  al.  (4),  has  been  modified  to  handle  the 
multiple-pulse  cases  of  interest  here  and  to  run  on  a  Digital  Equipment 
Corporation  microVAX  II  computer. 

The  model  is  divided  into  two  parts— one  that  computes  the  laser- 
induced  temperature  increases  in  the  eye  (for  spatial  coordinates  and  as  a 
function  of  time);  and  a  second  that  uses  the  predicted  temperature  In¬ 
creases  to  determine  the  degree  and  extent  of  irreversible  tissue  damage. 

The  mathematical  basis  for  the  temperature  predictions  is  the  standard  heat- 
conduction  equation  in  cylindrical  coordinates.  The  technique  used  is  to 
approximate  the  heat  conduction  equation  by  a  grid  of  finite  difference 
equations,  and  solve  these  at  successive  times  by  using  the  numerical 
solution  developed  by  Peaceman  and  Rachford  (16).  Since  the  resulting 
equations  are  linear,  temperature  rises  due  to  multiple  pulses  can  be 
predicted  by  adding  the  temperature  Increments  due  to  individual  pulses  with 
the  appropriate  temporal  offset. 
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The  time-temperature  history  predictions  are  then  fed  to  the  second 
part  of  the  model,  which  applies  the  Henriques  damage  criteria  to  determine 
the  extent  of  Irreversible  tissue  damage,  if  any. 

Appropriate  values  for  reflectivities,  transmissivities,  and  absorp- 
tivlties  for  each  ocular  component--as  well  as  geometric  data  and  indices  of 
refraction  for  the  rhesus  monkey  and  human  eyes— have  been  tabulated  (4,17) 
and  are  used  for  the  current  model  calculations. 


RESULTS 

The  ophthalmoscopic  lesion  tfireshold  was  first  determined  for  single 
100-ms,  647-nm  pulses  from  the  krypton  1  .or.  After  this  baseline  determi¬ 
nation,  thresholds  were  determined  for  two  identical  100-ms  pulses  with 
varying  times  between  exposures.  The  threshold  for  two  pulses  separated  by 
0.01  s  represented  the  shortest  Inter-pulse  interval  which  co.ild  be 
accurately  obtained  with  the  mechanical  shutters  used.  Other  inter-pulse 
intervals  were  0.06  s,  0.6  s,  6.0  s,  60  s,  and  600  s  (10  min).  The  ED50 
values  and  confidence  limits  for  this  series  of  experiments  are  listed  in 
Table  1.  (ED§o  used  to  designate  the  threshold  for  a  train  of  n  pulses 
in  terms  of  total  energy  delivered  by  Llit  puise  train.) 


TABLE  1.  RETINAL  THRESHOLDS  FOR  TWO  100-ms  PULSES  OF 
647-nm  KRYPTON  LASER  RADIATION(a) 


Pulse  separation 
.  .  (s) 

95%  confidence 
intervals  (mJ) 

o(c) 

2.6(c) 

2.3  -  2.8 

0.01 

5.0 

3.6  -  6.4 

0.06 

4.3 

...9  -  5.2 

0.6 

4.6 

4.2  -  5.1 

6.0 

4.5 

3.0  -  6.7 

an 

\j\j 

4.1 

'3  7  A  A 

w/  •  /  -  “T  .  ^ 

600 

4.4 

4.0  -  5,0 

_ 

4.5(d) 

4.2  -  4.8(d 

(a)  200-/im  retinal  image  diameter 

(b)  Total  energy,  two  pulses  . 

(c)  Single-pulse  threshold,  ED50 

(d)  Calculated  by  pooling  all  two-pulse  threshold  data. 
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A  seccnd  series  of  experiments  was  conducted  leaving  the  interpulse 
interval  fixed  at  0.06  s  (6.25  Hz  PRF)  and  varying  the  number  of  pulses. 
The  results  are  presented  in  Table  2. 


TABLE  2.  RETINAL  THRESHOLDS  FOR  100-ms  PULSE  TRAINS  OF 
647-nm  KRYPTON  LASER  RADIATION (a) 


Number  of  pulses  (n) 

(S3? 

95%  confidence 
intervals  (mJ) 

EDs^/pulse 

l(b) 

2.6(b) 

2.3  -  2.8 

2.6 

2 

4.3 

3.9  -  5.2 

2.2 

10 

20 

17  -  21 

2.0 

50 

87 

82  -  91 

1.7 

(a)  200-/im  retinal  image  diameter,  6.25  H^  FKF 

(b)  Single-pulse  threshold,  ED^q* 


All  of  the  foregoing  data  were  gener  ted  by  using  a  200-/tm  retinal 
image  diameter  of  the  laser  beam.  This  expanded  image  size  was  chosen  to 
insure  that  individual  pulses  in  a  multiple-pulse  sequence  could  accurately 
be  delivered  to  irradiate  the  chosen  patch  of  retinal  tissue.  In  order  to 
examine  the  dependence  of  multiple-pulse  additivity  on  retinal  image  size, 
the  threshold  determination  for  two  647-nm,  100-ms  pulses,  at  6.25  Hz  PRF, 
was  repeated  for  ~30-/4m  and  500-/im  retinal  image  diameters.  The  30-;im 
diameter  was  obtained  when  the  collimated  laser  beam  was  incident  at  the 
cornea. 


The  grid  of  macular  exposure  sites  was  adjusted  according  to  the 
retinal  spot  sire  to  avoid  overlapping  of  adjacent  grid  sitos.  Thus,  2x2 
grids  were  used  for  the  case  of  500-/im  retinal  image  diameter,  as  opposed  to 

200-pm  case,  and  4x4  grids  for  the  30-/im  spot  sizes, 
only,  4-6  additional  exposures  per  eye  were  delivered  to 
to  the  macula,  so  that  the  probit  calculations  could  be 
comparable  number  of  data  points  for  each  case.  No  significant 
was  found  between  the  macular  ED^p  threshold  and  the  paramacular 
ED50  for  either  the  single-pulse  or  the  double-pulse  exposures;  and  only  the 
combined  probit  results  (macular  plus  paramacular)  are  shown  in  Table  3. 


3x3  grids  for  the 
For  the  500-/im  case 
areas  just  adjacent 
based  on  a 


For  each  of  the  three  retinal  image  sizes,  the  single-pulse  threshold 
was  first  determined  for  comparison  to  the  corresponding  two-pulse  threshold 
(Table  3). 
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TABLE  3.  RETINAL  THRESHOLDS  EOR  lOO-ms  PULSES  OF  647-nm  KRYPTON  LASER 
RADIATION  AS  A  FUNCTION  OF  RETINAL  IMAGE  DIAMETER 


Retinal  image  diameter  Single-pulse  threshold  Two-pulse  threshold(®) 

{/im)  I  2 


“5? 

95%  confidence 
intervals  (mJ) 

95%  confidence 
intervals  (mJ) 

30 

1.0 

0.7  -  1.2 

2.3 

1.9  -  2.7 

200 

2.6 

2.3  -  2.8 

4.3 

3.9  -  5.2 

500 

6.0(b) 

:  0  -  6.6 

9.7(b) 

8.5  -  10.6 

(a)  6.25  Hz  PRF 

(b)  Calculated  by  pooling  macular  and  paramacular  exposure  data;  all  other 
thresholds  based  on  macular  data  only. 


DISCUSSION 


Experimental  Data 


The  data  of  Table  1  indicate  that  the  degree  of  additivity  for  two 
pulses  remains  approximately  constant  as  the  time  between  exposures  is 
varied  over  several  orders  of  magnitude  from  0.01  s  to  10  min.  When  the  two 
pulses  are  separated  by  >1  s,  thermal  equilibrium  is  certainly  reachieved 
before  the  second  pulse  is  delivered.  On  the  other  hand,  for  pulses 
separated  by  only  0.01  s  or  0.06  s,  the  irradiated  tissue  is  still  at  an 
elevated  temperature  induced  by  the  first  pulse  when  the  second  pulse  is 
delivered.  According  to  our  thermal  model  calculations  (discussed  in  the 
next  section),  the  temperature  at  or  near  the  center  of  the  irradiated  area, 
at  0.01  s  after  the  Initial  exposure,  has  relaxed  to  roughly  half  or  the 
peak  temperature  increment  achieved  at  the  end  of  that  first  pulse;  aiid,  at 
0.06  s,  the  temperature  has  relaxed  to  "'lOIX  of  the  peak  value.  Never¬ 
theless,  as  indicated  in  Table  1,  the  additivity  of  the  two  pulses  is 
Independent  of  whether  or  not  thermal  relaxation  is  completed  during  the 
interpulse  interval.  The  agreement  of  the  thermal  model  predictions  with 
this  observation  orovides  tentative  validation  of  the  numerical  olution 
technique  used  in  the  thermal  model  (4,16). 


This  result  appears  also  to  support  the  simple  empirical  model  that  the 
threshold  per  pulse  for  a  train  of  n  identical  pulses  is  related  to  the 
threshold  for  a  single  pulse  from  that  train  by  the  factor  n"!/^,  and  is 
independent  of  PRF  or  total  tr;in  length.  For  n  =  2,  the  n"l/^  factor 
indicates  that  the  threshold  per  pulse  is  5. 84  times  the  single-pulse 
threshold,  thus  equating  to  a  predicted  ED50  =  4.4  mj,as  compared  with  the 
observed  value  of  4.5  mJ  for  the  pooled  set  of  all  2-pulse  threshold  data 
from  Table  1. 
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A  second  conclusion  to  be  drawn  from  the  Table  1  data  Is  that  any 
tissue  repair  process  must  proceed  with  a  time  constant  »10  m1n.  We 
attempted  to  determine  the  additivity  for  two  pulses  separated  by  24  hr,  and 
found  a  threshold 'Indicating  little  contribution  from  a  repair  process  even 
over  this  time  frame.  However,  direct  comparison  of  these  data  with  the 
thresholds  listed  In  Table  1  Is  confounded  by  the  fact  that  the  single-pulse 
threshold  based  on  observation  of  a  minimum  visible  lesion  at  24-hr  post¬ 
exposure  Is  lower  than  that  based  on  the  1-hr  criteria  used  here  (18). 
Because  the  majority  of  published  retinal  threshold  data  (for  both  single¬ 
pulse  and  multiple-pulse  exposures)  Is  based  on  a  1-hr  criteria,  we  did  not 
pursue  further  the  question  of  pulses  separated  by  >1  hr.  Without  addi¬ 
tional  experimental  data,  we  suggest  that  the  retliial  repair  rate  is 
comparable  to  the  96-hr  time  constant  reported  by  Griess  and  Blankenstein 
after  exposures  to  458-nm  argon  laser  radiation  (14). 

When  the  data  of  Table  2  are  plotted  on  a  log-log  graph  of  threshold 
vs.  number  of  pulses  (Fig.  1),  a  linear  relationship  Is  found  with  a  slope 
of  0.91.  The  accepted  empirical  relationship  (7-9),  that  the  threshold  per 
pulse  for  a  train  of  n  pulses  of  pulsewldth  t  Is  equal  to  times  the 

single-pulse  threshold.  Is  equivalent  to; 

ED§o(t)  -  t.V4ED5o(t).  (1) 

which  predicts  a  slope  of  0.75.  Using  the  single-pulse  threshold  from 
Table  2  (ED^q  =  2.6  mJ),  the  predicted  values  of  ED§o  were  calculated  from 
Equation  1  and  plotted  on  Figure  1  for  comparison  with  the  experimental 
data.  (Further  comparisons  between  experimental  data,  empirical  model 
predictions,  and  thermal  model  predictions  are  included  In  the  next  report 
section.) 

The  data  of  Table  3  were  collected  to  assess  the  relative  additivity  of 
repeated  exposures  to  a  point  source  and  an  extended  source.  An  earlier 
study  on  multiple-pulse  thresholds  from  a  Q-swItched  neodym1um:yttrium- 
alumlnum-garnet  (Nd:VAG)  laser  (16-ns  pulsewldth,  1064-nm  wavelength)  found 
no  additivity  for  exposures  with  a  900-/tm  retinal  Image  diameter  (19);  but, 
for  collimated  beam  exposures  with  the  same  source,  the  multiple-pulse 
threshold  followed  the  empirical  model  (Eq.  1).  We  find  no  trend  indicating 
that  the  additivity  is  diminished  for  larger  retinal  spot  sizes.  The 
relative  additivity  is  shown  in  Figure  2,  where  the  two-pulse  thresholds, 
ED§o/2,  are  plotted  vs.  retinal  spot  size  for  comparison  to  the  cases  of 
zero  additivity  (I.e.,  ED50  =  2  ED50)  and  total  additivity  (ED50  =  ED50). 
Quantitative  assessment  of  additivity  for  the  collimated  beam  case  Is 
problematical  In  view  of  the  exposure  parameters  employed  in  this  experi¬ 
ment.  With  a  100-ms  pulsewldth  and  a  6.25  Hz  PRF,  the  duration  of  a  two- 
pulse  exposure  Is  260  ms.  Eye  drifting  of  sufficient  magnitude  to  cause 
significant  displacement  of  the  exposure  site  relative  to  the  30-/im  retinal 
Image  diameter  Is  certainly  possible.  If  not  probable,  within  this  time 
frame;  but  such  eye  drifting  would  be  undiscernible  by  the  pre-  vs.  post¬ 
exposure  fundus  camera  observations,  unless  the  drift  were  >50  /^m.  On  the 
other  hand,  with  a  200  /im  or  larger  retinal  image  dlameterT^a  drift  re¬ 
sulting  in  significant  displacement  of  the  exposure  sites  for  two  or  more 
pulses  would  be  easily  noticeable  by  fundus  camera  observation,  and  would 
serve  as  grounds  for  discarding  the  data  for  that  exposure. 
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We  believe  that  the  apparent  discrepancy  between  our  findings  and  those 
of  Griess  et  al .  (19)  may  be  related  to  the  damage  mechanism  prevalent  with 
each  set  of  exposures.  The  exposure  parameters  for  the  current  work  were 
chosen  to  minimize  the  possibility  of  Interference  from  anything  other  than 
a  thermal  damage  mechanism.  The  Q-switched  pulses  used  by  Griess  et  al , 
almost  certainly  infringe  upon  the  regime  where  photo-acoustic  or  mechanical 
shock  mechanisms  compete  with  the  thermal  mechanism;  and  the  shock  mechanism 
may  show  a  greater  dominance  when  larger  areas  of  tissue  are  irradiated,  if 
so,  then  the  thresholds  reported  by  Griess  et  al .  for  900-;im  retinal  spot 
sizes  cculd  result  from  a  peak-power  dependent  mechanism,  and  would  not  show 
additivity  of  repeated  sub-threshold  exposures.  For  the  minimal  retinal 
Image  sizes,  however,  the  thermal  mechanism  could  still  be  competitive,  if 
not  dominant,  thus  explaining  why  the  multiple-pulse  exposures  of  Griess  et 
al.  show  additivity  comparable  to  that  ...  c.t’i,rr  studies  (19). 


Thermal  Model ---Empirl cal  Model— Experimental  Data  Comparisons 

The  discussion  of  the  thermal  model  predictions  of  multiple-pulse 
thresholds  must  begin  with  an  examination  of  the  behavior  of  the  model  in 
predicting  single-pulse  thresholds.  In  fact,  the  thermal  model  has  been 
used  with  great  success  in  predicting  i  .'llnai  tii.csholds  for  single  laser 
exposures  for  wide  ranges  of  laser  beam  parameters,  Including  pulsewidths 
from  10"^  to  10^  s  (4,12).  We  choose  to  examine,  first,  the  variation  of 
threshold  vs.  pulsewidth,  according  to  model  predictions  compared  with  the 
dependence  derived  from  experimental  threshold  data  and  used  in  laser  safety 
standards  (5-7).  For  visible  and  near-infrared  wavelengths,  the  retinal 
damage  threshold  is  found  to  vary  as  t^/^,  where  t  is  pulsewidth  This 
empirical  relationship  holds  over  a  nominal  pulsewidth  range  from  T.8  x 
10“5  s  to  10  s  (or  longer,  depending  on  wavelength).  For  pulsewidths  much 
less  than  the  thermal  relaxation  time  of  the  exposed  tissue,  threshold 
should  be  independent  of  pulsewidth,  assuming  a  thermal  damage  mechanism. 


The  threshold  vs.  pulsewidth  behavior  has  been  plotted  in  Figure  3. 

The  solid  line  is  the  MPE  for  visible  wavelengths,  according  to  the  ANSI  and 
AFOSH  laser  safety  standards  (5-7).  As  already  Indicated,  the  slope  of  the 
curve  for  t  >  1.8  x  10”5  s  is  equal  to  0.75.  Thermal  model  predictions  are 
plotted  on  tfie  curve  for  two  cases,  both  of  which  assume  a  647-nm,  2-nini- 
diameter  beam  incident  at  the  cornea  with  a  Gaussian  profile.  In  the  first 
1-mrad  beam  divergence  is  specified,  and  the  model  uses  a  point- 
>1  tn  calculate  the  geometrical  distribution  of  light  on 
distribution  approximates  a  Gaussian  profile  with 
fim.)  In  the  second  case,  a  200 -/iin  retinal  image 
model's  point-spread  subroutine  is  bypassed),  to 
chosen  retinal  image  size. 


case,  a 
spread  functic 
retina.  (This 
diameter  ot  30 
specified  (the 
experimentally 
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a  l/e2 
diameter  i 
match  the 


For  both  the  collimated  and  expanded  beams,  the  thresholds  predicted  by 
the  thermal  model  approximate  the  MPE  curve  and  follow  the  appropriate 
linear  behavior  for  pulsewidths  either  much  less  than  or  much  greater  than 
the  thermal  relaxation  time.  For  long  pulsewidths,  however,  the  slope  for 
both  the  collimated  and  the  expanded  beam  curves  approaches  "'0.90,  compared 
to  the  value  of  0.75  defined  by  the  safety  standards. 


Now,  having  examined  the  behavior  of  threshold  vs.  pulsewidth  for 
single-pulse  exposures,  we  consider  the  multiple-pulse  case.  The  empirical 
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Figure  3.  Log-log  plot  of  predicted  retinal  threshold  vs.  pulsewidth.  Model  predictions  are  shown  for 
both  a  collimated  beam  (minimal  or  30-ym  retinal  image  diameter)  and  a  diverging  beam  (200-ym 
retinal  image  diameter)  incident  at  the  cornea.  The  solid  curve  shows  the  MPE  levels 
according  to  prevailing  laser  safety  standards  (5,7). 


model  expressed  by  Equation  1  provides  a  reasonable  fit  to  published  ex¬ 
perimental  multiple-pulse  threshold  data  (8,9),  and  has  been  adopted  for  the 
treatment  of  multiple-pulse  exposures  in  the  latest  ANSI  laser  safety 
standard  (7).  Since  the  single-pulse  threshold  varies  as  (t)3/4,  Equation  1 
is  consistent  with  the  case  of  n  contiguous  pulses  of  pulsewidth,  t.  That 
is,  since: 

ED5o(t)  =  c(t)3M,  (2) 

where  ED^o(t)  is  the  single-pulse  threshold  for  a  pulsewidth,  t,  and  c  is  a 
constant--then: 

EDsoCnt)  =  c(nt)3/4 

=  n3/4ED5o(t).  (3) 

Since  the  right-hand  sides  of  Equations  1  and  3  are  equal,  the  formulation 
of  the  ANSI  standard  based  on  Equation  1  is  equivalent  to  stating  that  the 
threshold  for  a  train  of  n  pulses  of  pulsewidth  t  is  equal  to  that  for  a 
single  pulse  of  duration  nt,  independent  of  PRF  or  pulse  train  length. 

In  Table  4,  this  safety  standard  formulation  (I.e.,  the  empirical  ri3/4 
model  for  multiple-pulse  thresholds)  is  tested  against  the  experimental  data 
collected  in  this  study  and  against  the  thermal  model  predictions.  In  this 
case,  the  thermal  model  predictions  provide  a  significantly  closer  fit  to 
the  experimental  daua  than  those  obtained  from  the  empirical  model.  As 
shown  In  Figure  1,  when  the  experimental  data  were  entered  on  a  log-log  plot 
of  threshold  vs.  number  of  pulses,  an  excellent  fit  to  the  data  was  found  by 
a  line  with  a  slope  of  0.91.  The  thermal  model  predictions  also  imply  a 
straight  line  with  a  slope  of  0.91,  whereas  the  empirical  model,  by  defini¬ 
tion,  yields  a  slope  c ;  0.75. 


TABLE  4.  RELATIVE  RETINAL  THRESHOLDS  FOR  TRAINS  OF  100-ms, 
647-nm  KRYPTON  LASER  PULSES («) 


n 


Experimental  data(^) 

EDgo/ED^O 


Thermal  model  predictions  Empirical  model 
ED§o/t050  ED§o/ED^50  =  n3/4 


1 

1 

1 

1 

2 

1.7 

1,8 

1.7 

5 

- 

4.2 

3.3 

10 

7.6 

7.9 

5.6 

50 

34 

35 

19 

IS! 

200-/im  retinal  Image  diameLer, 
From  Table  2. 

6.25  Hz  PRF 
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To  exercise  further  both  the  thermal  model  and  the  tripartite  compari¬ 
sons  of  the  preceding  paragraph,  we  have  examined  the  data  of  Ham  et  al .  for 
trains  of  up  to  one  million  pulses  from  a  laser  scanning  device  (20).  The 
wavelength  is  again  647  nm;  the  pulsewldth  is  40  fts;  and  PRF  varies  between 
four  values:  100,  200,  400,  and  1600  Hz.  The  log-log  plot  of  threshold  vs. 
number  of  pulses  for  Ham's  data  is  shown  in  Figure  4.  A  least  squares  fit 
to  the  data  yields  a  slope  of  0.70,  whereas  the  thermal  model  predictions 
which  are  also  plotted  on  the  figure  are  F|t  by  a  line  with  a  slope  of  0.92. 

Summarizing  the  results  of  the  two  preceding  paragraphs,  we  find  that 
the  thermal  model  predictions  yielded  a  much  closer  fit  to  the  experimental 
data  generated  in  this  study  than  did  the  empirical  model;  but,  when 
compared  with  the  data  of  Ham  et  al .  (20),  the  relative  goodness  of  fit  of 
the  two  models  was  almost  exactly  reverc'^d.  In  both  cases,  the  goodness  of 
fit  of  either  model  seems  acceptable  relative  to  the  experimental  vari¬ 
ability  of  the  threshold  data;  and,  at  least  for  the  ranges  of  exposure 
parameters  examined  here  (where  a  thermal  damage  mechanism  is  thought  to 
operate),  we  can  feel  quite  comfortable  In  predicting  multiple-pulse  thresh¬ 
olds  when  the  corresponding  single-pulse  threshold  is  known. 

The  following  salient  points  are  subject  to  the  condition  of  a  thermal 
damage  mechanism  being  operative: 

1)  Additivity  of  multiple  pulses  Is  observed  even  for  retinal  image 
sizes  as  large  as  500  ftm. 

2)  The  repair  or  recovery  of  retinal  tissue  following  sub-threshold 
laser  exposures  is  slow,  with  a  time  constant  of  the  order  of  several  days. 
Therefore,  the  additivity  of  multiple  pulses  over  any  time  fr^.me  up  to,  for 
example,  24  hr,  is  not  affected  by  this  repair  process.  Only  when  the 
effects  of  exposures  repeated  on  a  daily  basis  are  of  interest  do  we  need 
consider  modifying  the  predictive  models  to  account  for  tissue  repair 
processes. 

3)  Now  that  the  validity  of  the  thermal  model  has  been  established  for 
multiple-pulse  threshold  predictions,  the  mociol  can  be  e.'.ercised  over 
broader  ranges  of  exposure  parameters  for  coi' pari  son  to  the  empirical  n^/^ 
model  and  to  other  available  experimental  dat..  Such  calculations  could  be 
valuable  in  identifying  conditions  where  the  simple  empirical  predictor  for 
multiple-pulse  thresholds  might  produce  poor  results  and  compromise  safety 
considerations.  For  example,  preliminary  thermal  model  calculations  indi¬ 
cate  a  widening  divergence  with  the  empirical  model  as  PRF  appreciably 
exceeds  1  kHz  (i.e.,  as  interpulse  intervals  become  short  compared  with  the 
thermal  relaxation  rate).  The  high  PRF  regime  is  of  significance  to  mili¬ 
tary  operations,  since  a  number  of  high  PRF  laser  systems  are  already  in 
use;  so  further  thermal  model  calculations  are  in  progress. 
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Figure  4.  Log-log  plot  of  threshold  vs.  number  of  pulses  for  40-vS,  647-nm  krypton  laser  radiation  at 
high  PRFs  (100-1600  Hz)  (20).  The  best  fit  straight  line  to  the  experimental  data  (solid 
line)  has  a  slope  of  0.70.  Thermal  model  predictions  for  the  corresponding  exposure 
parameters  are  also  plotted  (a),  and  fit  by  a  line  with  a  slope  of  0.92  (dashed  line). 
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